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ON  THE  HYPOTHESIS  OF  A  QUASI- STATIONARY  SURFACE  LAYER  | 

IN  THE  THEORY  OF  THE  TURBULENT  DIFFUSION  OF  A  HEAVY 

AEROSOL  IN  THE  BOUNDARY  LAYER  OF  THE  ATMOSPHERE  | 

r  •  '  i 

4  Following  is  the  translation  of  "O  gipoteze  kvaai-  : 

etatsionamogo  prizeamogo  sloya  v  teorii  turbulentnoy  I 

diffusii  tyazhelogo  aerozolya  v  pograniohnom  sloye  ; 

atmosfery"  (English  version  above)  by  I»L.  Karol'  in  : 

Izv,  Akad.  Nauk  SSSR,  ser.  geofiz,  (Trans.  AS  USSR,  geo- 
phys,  ser.).  No.  5*  I960,  pages  753-761. 7 

1.  In  papers  [\J  and  /" 2  J  on  the  semi-empirical  theory  { 
of  turbulence,  we  examined  the  non- stationary  process  of  vertical, 
turbulent  diffusion  of  a  heavy  homogeneous  aerosol  from  an  in-  ; 
atantaneous  point  source,  placed  high  above  the  surface  layer  of  j 
the  atmosphere.  Since  the  exact  solution  of  this  problem  is 
tedious  and  of  little  use  for  analysis  and  numerical  calculation,! 
a  simple  approximate  solution  was  described  in  /~l„7t  in  which  j 
the  variation  of  the  coefficient  of  vertical  turbulent  diffusion  ! 
in  the  surface  layer  was  taken  into  account  in  the  boundary 
condition  for  the  problem  at  ground  level.  { 

•  ;  I 

*  ! 

The  exact  solution  of  the  problem  can  be  simplified  in 
another  way.  In  papers  /~5.7»  /  4 J  and  /  5_7  on  the  theory  of 
;  temperature  distribution  and  the  thermal  transformation  of  air 

•  masses  in  the  boundary  layer  of  the  atmosphere,  it  was  shown 
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I  that  the  temperature  ohange  in  the  surface  layer  can  be  consider-; 
ed  as  quasl-stationary  from  some  time  sufficiently  long  after  ^ 
the  beginning  of  the  process,  i.e#  we  can  ignore  the  partial 
derivatives  with  respect  to  time  in  the  corresponding  different-  j 
ial  equations.  : 

i 

4 

| 

In  this  paper  we  discuss  the  possibility  of  applying  the  j 
hypothesis  of  a  quasi-stationaxy  surface  layer  to  the  problem 
of  the  turbulent  diffusion  of  a  heavy  aerosol*  and  the  connect-  | 
ion  between  the  approximate  solution  obtained  in  this  way  and  j 
the  approximate  solution  found  in  /~1 J*  j 

In  flj  and  /~2  _7,  the  dimensionless  concentration  ?  (?.  S) . 
of  the  aerosol  is  defined  as  that  solution  of  the  equation  j 


ft  y  /r\  • 

ez  W  dZ  * 


ft;  CO 
li;  t>i, 


(1*1) 


which  satisfies  the  initial  and  boundary  conditions 

ff<P.t)«*»(C-it) 


q-*  0,  t-v  +  ee;  [K  (Z)  dqf  dZ 


“  0f  0  “  ft 


(1.2) 

(1*3) 


We  ahull  use  the  notation  of  /”l„7  &n&  Z~2„7s  t»  t.  to.  X* 
v=»«rpare»  respectively,  dimensionless  time,  the  vertical 
coordinate,  the  height  of  the  ’unevenness  layer’,  the  height  of 
the  source,  the  velocity  of  the  gravitational  subsidence  of  the 
aerosol,  and  a  parameter  characterising  the  reciprocal  action 
of  the  aerosol  and  the  surface  of  the  earth  ( $  *  0  -  reflection, 
»  o©  -  absorption  of  the  aerosol  by  the  earth).  The  height 
of  the  surface  layer  is  taken  to  be  equal  to  unity  and  to 


be  greater  than  X  >1*. 


Taking  the  diffusion  in  the  surface  la; 
stationary,  we  replace  (l.l)  by  the  equation 


layer  to  be  quasi- 


9  rm  .  ,  ^1” 
as K as; 


*=  o,  c®  <  c  <  i; 


,  «> i.  (I.la) 


In  this  ease  the  concentration  *%((%  0  of  the  aerosol  is 
found  by  solving  the  problem  (l.la)  -  (1*3)  in  exactly  the  same 
w&y  as  the  problem  (l.l)  -  (1.3)  was  solved  in  £\jj  and  £ 2j<, 
in  particular  by  using  a  one-sided  Laplace  transformation  for 
the  time  X  and  by  subsequent  asymptotic  representation  of  the 
Mellin  rotation  integral  for  large  valuee  of  X  and  with;  Xu  —  X“~  1- 
In  this  way  we  obtain  the  following  asymptotic  representation 
for  C)  fori  i <  C <x  * 


ylS~(Xt~VT)V4T  r  /  V2  \  1  i 

(t,  c)*“  __  xi)  t?  X>  ”  **'  <?  » (1.4) 


Irfle  K 


v/2-±  V +  Po\  1\~  (x2£2  ~  v*>/4; 


!  (hi  -f  w) 6C«  —  A|  $  XW1  —  0  in  —  8  •••} 

iy«  •  *  > 0;  <*■  *■  »  -  tMhIH 

.  i  v  **  0.  ! (1*5)  I 

The  solution  ^(xX)ot  the  problem  (l.l)  -  (1,3)  in  the. 
layer  i<£<X  also  has  the  asymptotic  representation  (1.4),  but 
the  function  ’ijy ip  ,  £c ,B)  i a  then  a  complicated  expression  (see 
£zj)  involving  the  modified  Bessel  functions  7v  (2  V/^), ,  Kv  (2  Y'Po)t 


*  In  £ 2,J  we  also  discussed  the  case  when  there  is  a  steady 
change  in  the  coefficient  of  vertical  turbulent  diffusion  in  the 
surface  layer,  IX  (£)*=£*«  t<i,  0<e<  2,  for  a'weightless 

compound 1  ( a>  s  Q  ) . 
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I  Let  us  consider  the  oasev»0#  If  we  take  Xj  t  to  he  | 

■bounded  and  let oo,  then  po  ~  Xj/4-r*  ->  0,  —  i  Vpo~*  0.  •  ®hen, 

:  by  using  the  expression  for  Bessel  functions  for  small  arguments  ,j 
it  is  easy  to  see  that  in  this  case  \  (f>9>  becomes  CfeX0'$] 
i.e,  the  principal  term  in  the  asymptotic  representation  of 
<f(x,  O  ,  the  solution  of  the  ’exact'  problem  (1.1)  -  (l»5)» 
is,  for  large  values  of  *C  ,  equal  to  the  principal  term  of 

the  solution  of  (l.la)  -  (1*3),  the  problem  with  a  quasi¬ 
stationary  surface  layer*, 

'  ' 

We  shall  now  take  V  >0,  and,  for  simplicity,  we  put  ^  0** 
Then  the  expression  will  take  the  form  (see  /~2J7)* 


'MPo.  o, 


Wv+i  (2  YPo)~^K  (2  VPo) 


V Pa? *-j-i  (2  Vfi)  —  (2  V Pt>) 


(1.6) 


— ya/ 4=^0,  k. 


v/2 


As  r-*<»  and  for  bounded  Tt ,  ,  since 
in  this  case  (1,6)  does  not  become  equal  to  Yvl  (pt,  0,  0)*ski/k9 


for  (x,  ,  the  solution  of  (l.la)  -  (1.3),  although 

as  T-*-oo  the  limits  of  these  expressions  are  equal, lim (p*,  0,  0)  = 
lim T¥j (p0,  0,  6)  =  1  .  It  is  easy  to  see  that  this  also  happens  when 
>•  0  and  for  any  values  of  the  parameter  yS» 


\ 


*  .There  ia  a  similar  result  for  a  'weightless1  compound  in  the 
ease  JT(C)  ■* C*,j0<t  <  2  in  the  surfaoe  layer. 

**The  possibility  of  ignoring  the  height  of  the  'unevenness 
layer'  is  examined  in  LlJ  and  /~2_7»  where  it  is  established 
that  if  this  is  done  the  solution  does  not  depend  on  the 
parameter  , 
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Thus,  when  we  are  concerned  with  the  diffusion  of  a 
heavy  compound  (  V  >  0)  the  principal  term  in  the  asymptotic 

representation  of  the  solution  <gCi’ (t,  O  of  the  *  exact* 
problem  (l.l)  -  (1.5)  does  not  become  equal  to  that  of  <ft 
the  solution  of  the  problem  with  a  quasi-stationary  surface 
layer,  for  large  values  of  X  ,  although,  as  ,  these 

solutions  have  the  game  limits  equal  either  to  aero  or  to  the 
stationary  distribution  of  concentration  in  height  when  the 
aerosol  is  reflected  by  the  earth's  surface. 

Thus,  the  simplifying  assumption  for  the  problem  (l.l)  - 
(1,5)  that  the  change  of  concentration  of  the  aerosol  in  the 
surface  layer  is  quasi-stationary  can  be  made  for  large  values 
of  X  only  when  V  »  0,  i.e,  only  for  a  ’weightless’  aerosol*. 

2.  To  calculate  the  time  taken  to  set  up  quasi-station- 
ary  state  in  the  surface  layer  when  V  «  0  we  define  the 
relative  error  B  of  substituting  ftJCx.  £)  by 
putting 

E{X,  C)  =  i-?<2)(T,  o/gfc*.  5)  (2.1) 


and  restricting  ourselves  to  the  case  £«  =0.  Putting  the 
principal  values  of  and  given  in  (1.4)  in^C/k*  ^ 

defined  by  (1.6)  with  V  0  and  'F„i  (/>„.  0.  £)~  — i,  af  ter 
simplification  we  have,  on  the  boundary  of  the  surface  layer 


E(x,  1)  ~/„(x)l[h  (r) +  /„(*)); 

The  requirement  \g\  <  €  gives  an  inequality  for  X t 

(2,2) 


*  f  57  -  f  5j  concern  particularly  the  turbulent  diffusion 
of  ’T/eightless  substances’  -  heat  and  humidity. 


Prom  tables  of  Bessel  functions  (  H'&Jj  )  find  that 
for  £*  0.1,  x<0.2.  Passing  to  dimensional  quantities  (see 
£lj  a n&  £  2 J)  w®  obtain  an  inequality  for  the  time  needed 
to  ®et  up  a  quasi- stationer;;/  state  in  the  surface  layer  to  an 
accuracy  of  up  to  10%?  t  >  tQ  (where  tQ  «  500,  5000  second®, 

the  height  of  the  souro®  being  H-'IQ,  100,  1000  m  respectively, 
and  where  Kg /L  Im./eecs  K£  being  the  average  magnitude  of  the 
coefficient  of  vertical  turbulent  diffusion  outside  the 
surface  layer?  and  L  is  the  height  Of  the  surface  layer)* 

Prom  this  it  follows  that  in  many  problems  concerning 
the  turbulent  diffusion  of  a  'weightless*  substance  in  the 
boundary  layer  of  the  atmosphere  the  process  in  the  surface 
layer  can  in  practice  be  thought  of  as  quasi- stationary  from 
the  beginning.  In  this  ease  we  can  suggest  the  following 
simple  method  for  solving  similar  double-layer  problems,  taking 
an  arbitrary  law  of  change  for  the  coefficient  K(£)  in  the 
surface  layer* 

In  the  layer  0  4  £  <  1  the  solution  of  the  ordinary 

differential  equation.  (1* la)  for  V  »  0  ,  which  satisfies  the 
second  condition  of  (1*5) »  i®  of  the  form 


X« 


mct, 

'  dl 


(2.5) 


where  is  a  parameter  independent  of  £  «  Eliminating  G^» 
from  these  equations  we  obtain  the  relations 

!  [Hl)  AU)1  -0-  1  lA  (2^) 

;  bf-^  L,-0'.  T-T+ims- 


«  on  the  boundary  of  the  surface  layer  (£  «  1,  K  (l)  »  l)« 

As  e.  consequence  of  the  equations  *=  fig*®  fdt  |c_ 

9|D(t  1)  (2.4)  can  be  used  as  the  boundary 
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condition  at  the  height  K,  ®  1  to  determine  e^f  (Y,  £)  »  i»®» 
to  solve  equation  (2# la)  in  the  half-plane  £  3*  1  with  the 
initial  condition  (1.2)  and  the  first  condition  of  (1*5 ) *  This 
problem  is  easily  solved^  since  for  £  ^  1  elation  (l.la)  ha® 
constant  coefficients*  Having  found  ve  find  the  J 

of  (2*3)  by  using  the  equation  «  dq^f  dt  »  an<*  f^cm 

(2.3) ' we  can  find  the  quasi- stationary  solution  (t£)  in 

the  surface  layer*  We  observe  that  the  parameter  of  (2«4) 
is  given  by  v  . 

i  =  (2.5) 

where  3)  »  ^  J  -  is  the  experimentally  determined 

total  turbulent  conduct ibility  of  the  surface  layer  (  £lj)* 

3.  In  flj  we  described  a  method  for  reducing  the 
problem  (l.l)  -  (1.3)  approximately  to  a  boundary  problem  for 
an  equation  with  constant  coefficients  (equation  (1*1)  with 
K  (^)  »  l),  by  transferring  the  effect  of  a  ohange  in  K  (£) 
in  the  layer  approximately  into  the  boxin&ary  condition 

(1.3) ' at  the  height  .  {*  Condition  (1.3)  then  takes  the 
form 

+  =  to).  (3.1) 


where 


1  /  V  ■+  (i /?  -  l/v)(CJ  -  i-CoKi-v)”1 

i  +  (ift  —  i)  Co  (i  —  In  So) 

1  /  {S  — -  In  Co  4*  £o  —  1 


v»  0;  1 


(3*2) 


The  initial  condition  (1.2)  and  the  condition  as  £-*-+« 

*  In  flj  the  value  of  B  is  given  only  for  $  »  06  ,  V  >  0*; 

however,  the  method  described  there  can  be  used  to  find  B  for 
any  /0  0  and  for  the  case  V  =  0* 
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are  unchanged. 

let  us  find  the  connection  between  the  solution  y  ( T,  £  } 
of  this  problem  and  the  solution  £>  i  «*  1;  2)  of  the  problea 

with  a  quasi- stationary  surface  layer.  In  £  1 J  we  found  the 
principal  term  of  the  asymptotic  representation  of  » 
analogous  to  (l»4)i 


i  'h(~,  t) 


-  vT)-/4r  r 

-1 _ _ _  . V)  __  CUZ~£*) 


2.V- 


h±lz l*]Fi  +  o 


72  —  7  —  Co* 


-“jo-  V  i? 


(?)  '  (3.3) 


j  Putting  #  **  0,  let  us  form  the  expression  AoK)™. 

j^j .using  the  principal  term  of  q^  given  in  (3*5) »  As 
X*  / 1'™'®’ 0  we  obtain  in  the  limit 

j  .  A* 

lim  At>(r)  /i/  [1  -f- 1/  (I  — £ft)]  =:  [1  /P  —  In^o]1  =■**>,  (3.4) 


where  B  is  defined  by  formula  (2.4)  for  K  (£)  »  £  .  This 

means  that  for  large  values  of  time  "T  the  solution  q^  (  ‘t'.  4.  ) 

satiefi6>s  the  condition  (2,4)  and,  therefore,  a  quasi-station-  i 
ary  stage  state  is  set  up  in  the  surface  layer  £fc<  £’<  \  for  q^*  j 

We  can  verify  at  once  that  as  *t'~*  to  and  x£)  I 

tend  to  sero  for  1<£<X  end  v  «  0,  although  from  equations  j 

(3«2)  and  (2.4)  with  (£)  =  £  ,  ] 


iin»4t«w*  +  c--c»_1 

^oo<?{2)  1  /  ©  -h  C  —  1 


(3.5) 


Thus,  for  large  values  of  *C  ,  and  for  r>?  ,  the  solution  j 
and  the  solution  q.^  (  2  )  of  the  problem  with  a  quasi- stationary 
surface  layer  are  equivalent. 

Prom  this  and  the  fact  that,  as  we  saw  in  section  l, 
is  equal  to  q(*C,4)»  the  solution  of  the  *  exact*  problem 
(l.l)  -  (l,3),  when  V  *  0  and  as  Z~+oo ,  it  follows  that 


q2('C\  £  )  and  q(r,£)  are  equal  for  ^  >  l  and  for  large  values  of  I 

In  conclusion  we  note  that  when  V>0  the  principal  terms 
of  the  solutions  q2(f »?)  and  q^^(T»£)  defined  by  the  expressions 
(3,5)  and  (1.4),  (1.5)  respectively,  are  no  longer  equal  as  X  ~J*‘ 
co  with  >  i  ,  For,  in  the  eame  way  as  for  (5*5)*  it  is  not 
difficult  to  find  that 

..  ?*_  _ C-i^.d-  2/(2/>~v) _ 

*"«  cf 7  **  S  -  i  -h  2  (p  - 8^)  /  V  $  +  6Co)  • 

This  limit  is  not  equal  to  unity  for  any  1  if  B  has  the 
value  given  by  equation  (3.2).  This  corresponds  to  the  deduct-  j 
ion  in  section  1  that  for  V>0  and  a3  •  ?  .  the  solution  q(  »  )  .  j 
of  the  'exact*  problem  (l.l)  -  (1.3)  does  not  become  equal  to  ] 

the  solution  (XtX>)  for  the  quasi- stationary  surface  layer,  j 

and  so  the  solution  q2(*r»^)»  nearly  equal  to  t  is  j 

different  from  q^2)(-r,£)  for  V>0and  large  values  of  ~C  ,  ■  j 

Thus,  our  simplifying  assumption  concerning  a  quasi-  i 

stationary  state  in  the  surface  layer  is  legitimate  only  for  j 
the  problem  of  the  turbulent  diffusion  of  a  'weightless1  j 

compound  (aerosol).  The  solution  of  the  problem  of  turbulent  j 

diffusion  of  a  heavy  aerosol  (  V  >0)  in  the  boundary  layer  of  j 

the  atmosphere  can  be  simplified  by  the  method  explained,  in  j£  iji 
of  substituting  the  coefficient  of  K(g)  in  equation  (l.l)  j 

by  unity,  and  the  boundary  condition  (1.3)  by  conditions  (3.1)  j 
and  (3*2)  at  the  height  X>  «*  ^T0. 


*  In  £l.J  we  established  only  that  the  flow  &  <0^9 <*£  +  v9  .of 
the  function  q  (Ti£)  was  nearly  equal  to  that  of  qg^,  £  )  at 
the  height  £  *  fe,  * 
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